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'^ '. 1 Introduction 



The potential group method is applied to the n-dimensional Coulomb- 
Rosochatius potential, whose bound states and scattering states are worked 
out in detail. As far as scattering is concerned, the S'-matrix elements are 
computed by the method of intertwining operators and an integral repre- 
sentation is obtained for the scattering amplitude. It is shown that the 
maximal superintegrability of the system is due to the underlying poten- 
tial group and that the 2n — 1 integral of motions are related to Casimir 
operators of subgroups. 



The Coulomb system in arbitrary dimensions is among the well-known and best 

studied exactly solvable systems in quantum mechanics |T]-[T^. It describes the 

Jv>( , dynamics of charged particles under the influence of the 1/r potential. This is 

}_( ' what we conventionally mean by Coulomb system, even if the true n-dimensional 

5t ! Coulomb potential satisfying the Poisson equation goes like l/r"~^, for n > 3. 

It has been shown that the bound and scattering states of a Coulomb system 
in n dimensions can be associated with the rotation group 5*0 (n + 1) and the 
Lorentz group SO (n, 1) , respectively [T- "4". In other words, for a fixed energy, 
these groups appear as invariance groups of the system. These symmetries for 
n = 3 were discovered by Fock[T3] and Bargmann[TS]. Moreover, it was pointed 
out in Ref. ^ that the dynamical group that yields the positive and negative 
energy spectra of the Coulomb problem in n dimensions is SO (2, n -I- 1), while 
the dynamical group in the case n = 3 had been discussed by Barut[TH]. It is 
also worth noting that Coulomb and harmonic-oscillator systems are the only 
spherically symmetric superintegrable systems in arbitrary dimensions. 



We recall that in classical mechanics a closed system with n degrees of free- 
dom is completely integrable if it admits n integrals of motion (including the 
Hamiltonian ) that are independent and in involution , i.e. the Poisson brack- 
ets of any two integrals are zero. The system is called superintegrable if there 
exist q, 1 < q < n — 1, additional independent integrals of motion. The cases 
q = 1 and q = n — 1 correspond to minimal and maximal superintegrability, 
respectively. In quantum mechanics the definitions of complete integrability and 
superintegrability are same, but Poisson brackets are replaced by commutators 
(see Ref. [T7] for a general review) . 

In the present work we shall consider a quantum system characterized by a 
Hamiltonian of the form 

i—l * 

We shall call the system governed by the Hamiltonian given above the 
Coulomb- Rosochatius [18] system. The motion is confined to the region bounded 
by the singularity of H at hyper-planes Xi = 0, i = l,2,...,n , i.e. hyper- 
octant ( one of the 2" regions of Euclidean space £'"). Without loss of gen- 
erality this can be taken to be the non-negative hyper-octant ( i.e. where 
Xi > 0, i = 1,2, . . . ,n ). 

The main interest of the Coulomb-Rosochatius system consists in its maxi- 
mal superintegrability, recently proved for n = 3 in the classical casejlBj (see also 
[20]). The system had been proved long ago [21], [22] to be quasi-maximally su- 
perintegrable for n < 3, i.e. to admit 2n— 2 operators quadratic in the momenta. 
For n = 3, the fifth integral of motion quartic in the momenta was explicitly 
worked out in Refs[T9|-[2^, thus proving the maximal superintegrability of the 
classical Coulomb-Rosochatius system. Finally, the maximal superintegrability 
of the classical system ([T]) in n-dimensional spherical, hyperbolic and Euclidean 
spaces was proved in Ref. |23) : here again, one of the 2n — 1 functionally inde- 
pendent integrals of motion turns out to be quartic in the momenta, while the 
remaining ones are quadratic. 

The main purpose of the present work is the complete algebraic solution 
of bound states and scattering states of the quantum-mechanical Coulomb- 
Rosochatius system in n dimensions and the derivation of the 2n — 1 integrals 
of motion expected for maximal superintegrability. 

Before describing the method of solution of the quantum Coulomb-Rosochatius 
system in n dimensions , a few definitions are in order; a Lie group G is an in- 
variance group for a quantum-mechanical system with Hamiltonian H if the 
latter can be related to a suitable function of a Casimir operator, C, of the 
group, G 

H = f{G) . (2) 

An algebraic derivation of the energy spectrum is possible also when 

H = f{C)\n. (3) 



where ?^ is a subspace of the carrier space. In this case the group G describes 
the same energy states of a family of Hamiltonians H with different potential 
strength. This is why G is called potential group [53]. Such an approach was 
proposed by Ghirardi [3S] , who worked it out in detail for the Scarf potential 
pS] . It is similar to the approach of Olshanetsky and Perelomov [371 HB] ■, where 
quantum integrable systems are related to the radial part of the Laplace op- 
erator on homogeneous spaces, i.e. to the radial part of a second-order Casimir 
operator of a Lie group. Moreover, it has been shown in Ref. [15] that the S 
matrix can be associated with an intertwining operator, A, between two Weyl- 
equivalent representations, U^ and C/-*^, of G, i.e. two representations with the 
same Casimir eigenvalues. 

By definition, A satisfies the following equations 

AU^{g) = U^{g)A,'igeG (4) 

and 

AdU'^ (b) ^dU'^ib) , V & e (5) 

where dU^ and dU^ are the corresponding representations of the algebra, g, of 
G. Eqs. (|3||5]) have high restrictive power, determining the intertwining operator 
up to a constant. The S matrix coincides with the intertwining operator A 

S^A (6) 

if eq. (0) holds, and with the reduction of ^ to a proper Hilbert subspace % 

S = A\h (7) 

if eq. © holds. 

The plan of the paper is as follows: Section 2 will describe in full detail the 
potential group for the Coulomb-Rosochatius system, with a potential special- 
ized to 

VM-J + E^^^^. (8) 

with Subsection 2.1 dedicated to the derivation of bound states and Subsection 
2.2 to scattering states. Finally, Section 3 will be devoted to conclusions and 
perspectives. 

2 Potential group for the Coulomb-Rosochatius 
system 

It is well known that the group SO {N + 1) ( SO {N, 1) ) has a class of unitary 
irreducible representations (UIR's) characterized by a number j, j = 0, 1, 2, . . . 
(j — — ^^^ + ip, p > 0) in which the basis vectors are completely labelled by 
iV — 1 numbers only. This representation can be realized in the Hilbert space 
spanned by negative-energy (positive-energy) states corresponding to a fixed 



eigenvalue of the Coulomb Hamiltonian in N dimensions. According to this, 
one introduces the angular momentum operators L^ and the A'^-dimcnsional 
Runge-Lenz vector Ai 

Li] = C^T^j - OtTj , {i,j = 1, ■•■, N) (9) 

and 

JV 



^^ = -T^'^iLiji: 






(C • ^) TT, - Gvr^ - J^^TT, + ^, (* = 1, ..., N) (10) 



respectively, where C = (Ci. C2, • • • , Qn), tt = ( tti, 7r2, . . . , ttat) , ttj = -«^, , 

C^ = X]i=i C^' ■"'^ — X]i=i ^i^- These operators satisfy the following commuta- 
tion relations 

[Lij,Lki] ~ i {SikLji + SjiLik — SiiLjk — SjkLu) 

[Lij,Ak] = i {SikAj - SjkAi) (11) 

where h is the Coulomb Hamiltonian in N dimensions 

h^\^-^. (12) 

2 ^ ^ ' 

Since [X^ , /i] = \Ai , /i] = , we may restrict the above algebra to a subspace 
where h has a definite eigenvalue e and define a new set of operators as follows 

M,,o = -Mo,»-|2er^A, , (z = l,...,iV) 
M,, = Ly , (i,j-l,...,iV). (13) 

Then, as a result, we obtain the Lie algebra of SO {N + 1) (when e is negative) 
and the Lie algebra of SO {N, 1) (when e is positive) 

[M^^, M„x] - i (g^caM^x + 9uxM^„ - g^x^U^ - g^^M^x) , (14) 

where iJ.,^ = 0, 1, 2, ..., N and 

5^^- = (+,+,...,+,+) for SO{N + l) (15) 

5m>> = (-,+,...,+,+) for S0{N,1) 

The generators -M^^ act in the eigenspace "H of ft, equipped with the scalar 
product 

,h)= I <t^l{0^2{0dC Ce^"^ (16) 






where dC = dC,idC,2 ■ ■ ■ dC,N- A detailed discussion of the SO {N + 1) [SO {N, 1)) 
representation generated by operators (fT3|) is given in Refs.[T]-[4]. 

At this stage we note that, in general, one can define the generators of 
SO (N+l) ( SO {N, 1) ) (let us call them M^^), as follows 



M; 



1,0 



M,> 



= -Moa = |2er^ Ai/2 (^) o A, o X'^/^ (C) 



where A is some non-negative function of (. Now the generators A/^i/ act in the 
eigenspace T-L oi h — A^^^ {Q) oho A^^/^ {(,) equipped with the scalar product 



(01,02J 



<^t (C) -^2 (C) rfA^ (C) , C^R 



N 



(17) 



where d/i {(,) = A^^ (C) dC, is a quasi- invariant measure on R^ . The represen- 
tations acting in % and "H are, of course, unitarily equivalent. The unitary 
mapping W which realizes the equivalence is given by 



W 



^^0 = Ai/2(^)^ 



(18) 



Although these representations are equivalent from the mathematical view- 
point, they may be related to different physical problems. We shall prove that 
the bound states and scattering states of the quantum system ([5|) are related 
to the representation of 5*0(3^ -I- 1) and 5*0(3^, 1) , respectively, acting in the 
Hilbert space with scalar product 



(^1,^2)= f ri (0^2(0 dfiio 



(19) 



_R3" 



where d^ (C) ~ ]Xi=i {Cit-2 + Cli-i + Cli) d(. To this end, we choose N — 3n 
and A (C) — n"=i (C3(-2 + Cli-i +C3i)- Then we introduce the second-order 
Casimir operator 



-J on 



(20) 



and consider the following operator depending on C 

1' 






d 



^ d d d 

2 /-2 , a2 , a2 \ QiTTT + Q2^rr + Qs. — 



C3n-2 + Csri-l + Csn 



Csn- 



d_ 

3TI-2 



■50 



'dC2 

+ C3n- 



ac3 

d 



'TQ: 



3n-l 



+ C3n 



dC.'in 



27 



Decomposing i?^" into the direct sum of three-dimensional subspaces and in- 
troducing spherical coordinates in these subspaces 



C = {xiUi,X2U2,-..,X„Un) 



(21) 



where Uj = (sina; sin/3i,sinai cos/3i, cosa^) , Xi > 0, < a^ < tt, < /3i < 
27r, {i — 1,2, . . . ,n) , we have 



7 



[c+i'-^y 



(22) 



E 



92 



1 



1 



d_ . _d_ 

dxf xj \ sin Oi dai dai 



■ 2 
sm a, 



dP! 



27 



with X = {xi, X2, ■ ■ ■ ,Xn) ■ With this coordinate system the measure d/i (C) in for- 
mula (J19p becomes d/i (C) = rfa; Ilj^i sin aidaid/3i, where dx = dxidx2 ■ ■ ■ dxn- 

Let Hk, K = {ki,K2, . ■ ■ ,Kn',cri,a2, ■ ■ ■ jCTn), be a subspace of functions 
(j){() = * (x) nr=i ^Ki' (o^i: /^i) : '^i^h hxcd K, whcrc Y^.' (ofi, /3i) are spherical 
harmonics of degree k^. Thus, the operator ([^ restricted to this subspace 
becomes 



c+i'-^) 



n„ 



E 



dxf 



Ki (Kj -t- 1) 



27 



(23) 



where we have used the fact that 



d 



d 



-. "5— smaj-— 



Hence, the Hamiltonian 



sm a., 



dPt 



Y^l {a,, 13,) = -K, («, -t- 1) Y^; (a„ ft) 



2 



■^ Kj {Hj + 1) 



2x? 



(24) 



can be described in terms of the potential groups SO (Sn + 1) and SO (3n, 1) 
since 

H 



7 



3n-l\2 



c+{^) 



-Hi. 



Moreover, since the Casimir operators of the groups in the chain 

G D SO (37i) D SO (3?i -3)xSO (3) D 5(9 (3n -- 6) x S-O (3) x SO {3) D . . . 
D SO{3)xSOi3)x ...xSO{3)d SO{2)xSO{2)x ...xSOi2) (25) 

where G is either 5*0 (3n -I- 1) or SO (3n, 1), form a complete set of 3n com- 
muting operators ( including the Casimir operator, C, of G ), i.e. the opera- 
tors C, (7SO(3n-3p)^ C-SOCS), ^^^ (7^0(2)^^ ^j^g^.^ p = Q, 1, 2, . . . , n - 2, fc = 



1,2, ... ,n and 



/jSO{3n~3p) 



3n—3p 






/ ^ -7.2 



2 -^^-^ V (9a;, 9xi 



(26) 



1 d 



1 92 



■ sma 



sin aj daj daj sin aj dfij 



— 2 (n — p) (n — p — 1) , 



3(n-fe)+3 
(.50(3). ^1 ^ ^2 



■ij— 3(n — A;) + l 



19. d 

■ smafe- 



1 92 



sinafc duk "'"'"" 9afc sin^ afc S/S^ / ' 



and 



3(n-fc)+3 



(27) 
(28) 



nSO{2)k _ A V /if 2 - __ 
^ "2 ^ ^^^^^ dpl 

i,j=3(n-k)+2 '^ 

are mutually commuting operators, it follows from formula ([26l) that the n — 1 
operators 



Ip = 



r(^50(3«-3p) + 2 (n - p) (n - p - 1)1 



iVfr— _r— 1 V ^2 V Kj («J + 1) 

4=1 3 = 1 3 



2 '^—' \ ' dxj ^ dx. 

1,3=1 ^ ■> 



(29) 
, p = 0,l,2,...,n-2 



are integrals of motion. These integrals of motion are responsible for the sepa- 
rability of H in spherical coordinates. 

The remaining 71 — 1 integrals of motion can be related to Casimir operators 

ij=3q+l i,j=q+l \ J '/ yi=5+l J 



n ^ 

j=q+l 3 



Id. d 1 d^ 

■ smck,^:^ h 



sin a J duj ^ duj sin^ a^ d(3j 



2{n^q)in-q-l), 



C^ 



,3 , n 



fi,iy—0 



2e 



1 92 



Id. d 

■ sm a, 



p • x—yX ■ p — 2e 



Id. a 

■ smai- 



1 92 



sinai 9ai 9q;i sin ai 9/3i 



(n-l)(n-3) 



4a;? 



H 2 



2£ 



withal = (a; ■ p)pi-xip'^-i^pi+^, x = (a;i, X2, . . . , x„) p = { pi,p2, . . . ,Pn) 
and pj = —i-J^, of the groups in the chain 

G D 50 (3n ~3)xG' D SO (3n - 6) x 5*0 (3) x 50 (3) D . . . 

D 50 (3) X 50 (3) X ... X 50 (3) D 50 (2) x 50 (2) x . . . x 50 (2) 



where G' is cither SO (4) or SO (3, 1). Namely, 



J, = 



and 



^SO(Zn-M) + 2 (n - g) (n - g - 1)] 
1 ^ / _d__ _d_ 



■H^ 



'ij=9+l 






Kj (Kj + 1) 



, 9 = 1,2,.. 



(30) 

2-2 



^n- 



= -2e 



Ki (Ki + 1) 



■Hf 






+ Ki (ki + 1) (p • a;) ^ (x • p) 



(31) 

(n-l)(n-3) 
4^1 



are constants of motion, too. 

It is worth noting that, since Jn-\ is quartic in the momenta, the com- 
plete set of mutually commuting operators {i/, Jg, g = 1, 2, . . . , n — 1} does not 
specify a separable coordinate system. 

2.1 Bound states 

The bound-state spectrum is immediately obtained from the eigenvalue of the 
Casimir operator C of the potential group 5*0 (3n + 1), i.e. j (j + 3n— 1), in 
the form 



E- 



7 



Hj 



3n-l) 



(32) 



where j takes on integer values from ki + K2 + . . . + Kn upwards. 

The basis functions \j; IMK) for Hilbert space Ti. can be defined as the 
common set of eigenfunctions of the Casimir operators of the groups forming 
the chain (1^ 



j{3 + -in-l)\3;lMK) 

rup {nip + 3n — 3p — 2) \j; IMK) , p = 0, 1, . . . , n 

K,{Ki + l)\j;lMK), i = l,2,...,n 



C\j;lMK) 

CSO(3h\^.jMK) 
CSOi2hy.jMK) 



a'i\i]lMK), i = l,2,...,n 



where mg = I and M and K are the collective indexes {mi,m2, . . . ,mn-2) 
and K = (ki, K2, ..-,'*«; cti, cr2, ..., a„), respectively. It is also known that in 
correspondence with every chain of subgroups in SO{N) it is possible to define 
a polyspherical system on R^ (see Ch. 10 of Ref.[5D]). According to this, we 
introduce polyspherical coordinates on i?'^" by choosing Xi in (PT|) as 



Xi = r sm t/n-i sm 6'„_2 • • ■ sm U2 sm t/i 
X2 = rsin0„_i sin0„_2 ■ • ■ sin02 cos^i 

a;„_i = rsin6'„_iCOs6'„_2 
a;„ = rcos9„-i 



where r > 0, < 6*, < 7r/2, for i = 1, 2, . . . , n - 1. 
By constructfon 



(33) 



i=l 

where ip (x) is the bound-state wave function 

ip (x) = TZji (r) J^jAf (x), X ^ x/r 



(34) 



Here TZji (r) is the radial part of the wave function, while yiM (x) is the angular 

f ii 
the 3n-dimensional Coulomb wave function [5] as 



part of it. According to p^ TZji (r) is related to the radial part 7?.^°"' (r) of 



7^,,(r) = r"7^,^''«' (r) 
where L^ are Laguerre polynomials and 

4(n+l) 



= (27) 



-n/2 



J- + ^(3n~l) 



r(j-^ + i) 

2r {j + l + 3n- 1) 



(35) 



while yiM (x) are related to the (3n— l)-dimensional spherical harmonics Yimk ( C 
in the above polyspherical coordinates as 

n— 1 n 

YiMK (c) = J^iM (£) n (sin""^0„-fccos0„-fc)"'[]l7; (a„/30 , C = CA 

fc=l 4=1 

(36) 
In order to write Yimk ( C ) "^6 make use of a graphical method called the 

tree method (see Section 10.5 of Ref.[30])- According to this method polyspher- 
ical coordinate systems on the sphere of unit radius can be described by graphs 
called trees. The trees contain nodes and each node has two edges, which are 
distinguished as left or right edge. An angle is associated with each node and the 
sine (cosine) of this angle is associated with the corresponding left (right) edge. 
The free (upper) ends of tree are labelled left to right with Cartesian coordinates. 
Then the Cartesian coordinate Q is equal to the product of trigonometrical func- 
tions on the edges along the unique path connecting the lowest node with Qi- In 
Fig.l we associate a tree with a polyspherical coordinate system on S*"^"'^ given 




Figure 1: A tree of polyspherical coordinates for S 



3n-l 



by formulae (|2ip and (|33)). Hence, the (3n— l)-diniensional spherical harmonics 
YiMK ( C ) corresponding to this tree can be obtained by the rules described in 



Section 10.5.3 of Ref. j30j . As a result we have 

n-2 



yiM (x) 



(m,+ ='"~^' -l, K, + i) 



X n «-"■+"'' Or.., C0S-+1 0„_.P(r;::-:.-; 



)/2 



(cos26'„_i) 



i=l 



xsin''"+^6'iCos'""-i+^6'iP, 



(k„ + ^, k„-i + ^) 



>(",/3) 



, w„ (cos 26*1) , (37) 

where P^'^' are Jacobi polynomials and x is a normalization constant 

r (i (rrii-i + TTii + Ki + 3n — 3i + 1)) F (^ (mi_i — to.; — k^ + 2)) (2mi + 3n — 3i + 1) 



n 



F (i (mi_i + r7i,j - Ki)) V (i (mj_i - m, + k,, + 3)) 
F {\ (m„-2 + «■« + Kn-l + 4)) F [\ (TOn-2 - K^ - K„_i + 2)) (2m„_2 + 4) 



F (i (to„_2 + k„ - K„_i + 3)) F (i (m„_2 - Kn + K„_i + 3)) 



(38) 



2.2 Scattering states 

Once the group structure of the problem has been recognized, the associated 
S matrix can be computed by using matrices that intertwine Weyl-equivalent 
representations of 5*0 (3n, 1) in the bases corresponding to the reduction (l^5|) . 
We find it expedient to use, for this purpose, equation @. By realizing the 
principal series of SO (3n, 1) on suitable Hilbert spaces of appropriate functions, 
one can derive from ^ the functional relations satisfied by the kernel of the 



10 



intertwining operator, written in integral form and, consequently, the explicit 
representation of the matrix elements of the operator itself. 

It is known that the most degenerate principal series representations of 
SO{N,l) labelled with the quantum number j — —■^^5-'- + ip, with p > 0, 
can be realized on £2 (S"^"^) (see Section 9.2.1 of RefjSU]) 

f/,(5)/W = M'7(ry,) , veS^-' (39) 

where 



W, 



a 



E-1 , , ^ l^t=l 9ki Vi + 9kN 
9mV^ + 9NN, {Va)k = ^^Iv \ 



i=l 



'Ei=l9N;Vr+gNN 



The representations specified by labels j and 1 — N — j are Weyl-equivalent. 
The operator A defined by 

{Af)iTj) = JlCiv,v').fiv')dv' (40) 

intertwines representations j and 1 — A^ — j on condition that 

/C (^„ rQ = {u,r-'+^ [u^'f-'^^ JCirj, ry') . (41) 

The kernel, K , is uniquely determined by Eq. (|4ip up to a constant and is 
given by 

/C(7?,r;') = ><(l-r?.ry')'"'^^' ■ (42) 

with 

_,_^,.,rgi+M ,43) 

TT 2 r [—ip) 

Taking into account the fact that the spherical harmonics Yimk P6p form a basis 
in £2 (S''^""^), corresponding to the above reduction, we obtain the following 
integral representation for the matrix elements of A 



{j- I'M'K'l A \j; IMK) = j IC{rj, v')Y;^i'K' iv') Yimk (v) dildr^' . (44) 

Therefore 

(j; l'M'K'\ A \j- IMK) = Ai5u'5mm'5kk' , (45) 

where 

r(M + zp + 

According to this, we have 

Sip;p')^Y. ^ly^M (P) ytu if) ■ (47) 

IM 



11 



Thus, the scattering amplitude, / {p',p')^ is defined by 

n-l 

f {p;p') = H) (—) ' E (^' - 1) ^"^ (p) yiM ip') ■ (48) 

^ P ^ IM 



We can omit unity in the brackets of formula (1481) when p' ^ p , leaving 

/ ^ r / 3n- 

;a/ V 2 
Moreover, formulae ([55)) and the following expansion of the kernel 



/(p;p')=(-.)(y) " E r(IS!ip+!) y"'''"y»'<^''- <« 



(1-,, ,/)---■-., 2,,^ jTg^i: £(s^^',w(,')w(,) , 

^ V 2 ^ 'A^j (A/if '- \ 2 H'^ '') 

yield an integral representation of the scattering amplitude 

1 Oior / 3n-l I ■) 

fip;p') = ^^ i^^^A(0„_i,...,^i)A(0U,...,0l) (50) 

ip 2 i (, — JPj 

X / ••• / (1 -pip'jcosai -p2P2COSQ;2 p„pJ^cosa„) 2 *'' 

Jo Jo 

xPki (cosQfi) P„2 (cosQ!2) • • • Pk„ (cosa„) sinai sina2 • • • sina„dai(ia2 ■ ■ • dctn 
where 

n-l 



A(0„_i,...,0i)=n 



sin" ' 6'„_i cos 6'„_i 



1=1 



When Ki = 0, formula (l50l) simplifies to 
2 2*T ("~^ + ip) .^-^ 

/ {P;P') = -^T^ F7^— 1 E ^ (1 ~ eiPlP'l - e2P2P2 (-nPnP'nY 



-ip 



(51) 



with 



n< 



j=i 
It is worth noting that the amplitude (ISTj) does not reduce to the Coulomb 
amplitude fcoui in fi dimensions 4 

JCoul {P;P) = —;:^ J., . . (1-PlPl -P2P2 PnPn) ' 

ip^r- 1 {-ip) 

(52) 
when Ki is set equal to zero. The reason for this discrepancy lies in the fact that 
we solve the Schrodinger equation for Xi > 0, (« = 1, 2, . . . , n) with boundary 
conditions 

* (a;) = for a;^ = 0, (i = 1, . . . , n) 



12 



3 Conclusions and outlook 

The present work is the latest in a series of papers where the potential group 
approach and the method of intertwining operators (formulae ([2][7l)) have been 
applied to non-central extensions of the Coulomb potential: Ref.^Sl studied 
the bound states of the three-dimensional Coulomb potential plus a barrier 
term with the SO (5) potential group and Ref. J3^ the scattering states of sim- 
ilar potentials with SO (5, 1) potential group. A first simultaneous analysis of 
bound and scattering states of a three-dimensional Coulomb-Rosochatius po- 
tential of type ([U was performed in Ref [33] , where the use of potential groups 
SO (7) for bound states and SO (6, 1) for scattering states. Generally, in the 
n-dimensional case for Coulomb-Rosochatius potential one can choose poten- 
tial groups SO [qn -\- 1) and SO {qn, 1) with q = 2,3, . . . . Then the potential 
strengths f3i will be related to eigenvalues ki {ki + q — 2) , ki = 0,1,2, ..., oi 

Casimir operators of subgroups SO {q)^ as: (3i — (ki + ^^) —j.ln the present 
work, we have chosen q — 3. 

The method has proved useful also in deriving the full set of 2n—l constants 
of motion, related to Casimir invariants of subgroups appearing in the decom- 
position chains, thus proving the maximal superintegrability of the system. 

The potential group approach is quite general and, obviously, not limited 
to the orthogonal and pseudo-orthogonal symmetries underlying the Coulomb- 
Rosochatius Hamiltonian. Systems with different symmetries will be studied in 
future works. 
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